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abstract. Molecular analysis of parapoxvirus envelope genes was performed. Parapoxvirus DNA was detected in eight calves from eight 
farms in Iwate Prefecture, Japan, between April and September 2010. Seven of the detected viruses were identified as bovine papular 
stomatitis virus (BPSV) by sequencing, because their nucleotide identity was more than 96.8% similar compared with BPSV strain V660. 
Among them, two formed a subgroup, because their amplicons were digested with Xmn I (a marker for BPSV) and Hinc 11 and exhibited a 
T61C nucleotide substitution in the sequenced region. The remaining virus was pseudocowpox virus that had not been reported previously 
in Japan. Our results demonstrate the presence of a new BPSV variant in Japan with genetic variability in the envelope gene. 
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Parapoxvirus infections are widespread in ruminants, usu- 
ally inducing mild clinical symptoms around the mouth or 
on the teats of affected animals [3]. These infections need 
to be differentiated from other serious diseases affecting 
the oral cavity and skin, such as foot-and-mouth disease, as 
clinical symptoms are similar [1, 6]. 

Known parapoxviruses (family Poxviridae) include orf 
virus (ORFV), bovine papular stomatitis virus (BPSV), 
pseudocowpox virus (PCPV) and parapoxvirus of red deer 
in New Zealand (PVNZ) [12]. Parapoxviruses were formerly 
classified based on natural host range, clinical symptoms and 
serology [13]; however, these criteria do not always reflect 
the molecular characteristics of these viruses [8, 17]. 

Amplification by polymerase chain reaction (PCR), 
sequencing and restriction fragment length polymorphism 
(RFLP) have been used for molecular characterization and 
classification of parapoxviruses [7, 8]. These methods have 
made it practical to study the prevalence of parapoxvirus 
infections in animals and to differentiate infections by other 
pathogens. Our aim was to determine unique molecular fea- 
tures of parapoxviruses detected in calves between April and 
September 2010 in Iwate Prefecture, Japan. 

Affected calves (n=14) were 2- to 10-month-old and 
found on eight farms (farms 1-8). Calves presented with le- 
sions (papule, erosion or ulcer) in the mucosa of the gingiva, 
lip, palate or tongue and in rare cases, excess salivation. No 
lesions were observed on teats or in the interdigital cleft 
of limbs. Clinical signs also included mild pyrexia (39.6- 
39.9°C) in three calves. Swab samples (n=10) were collected 
from mucosal lesions of eight affected calves and from the 
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gingival mucosa of two calves without clinical symptoms 
that had been housed near calves with clinical signs. There 
were no obvious differences among the clinical features 
found in 14 affected calves. The eight farms were scattered 
throughout the prefecture, although two farms (farms 6 and 
7) were located in the same town. There was no direct con- 
tact among the cattle and/or people between each farm. 

PCR, RFLP, sequencing and phylogenetic analyses were 
performed on the swab samples [7, 8]. Viral DNA was 
extracted with a QIAampDNA Minikit (Qiagen, Tokyo, 
Japan), followed by PCR with primers (PPP-1 and PPP-4) 
[7] to amplify partial nucleotide sequences (594 bp) of the 
envelope gene. RFLPs were conducted with Drd I (marker 
for ORFV), Xmn I (marker for BPSV), PflM I (marker for 
PCPV) or Hinc II (marker for PVNZ). Following direct nu- 
cleotide sequencing, the deduced amino acid sequences were 
aligned, and phylogenetic analysis was performed using the 
Maximum Likelihood Method. Phylogenetic trees were con- 
structed with the assistance of MEGA5 software. Nucleotide 
and amino acid sequences obtained were compared with 
each other with BPSV (V660; GenBank Accession No. 
AB044793), PCPV (VR634; AB044792) and PVNZ (DPV; 
AB044794) reference strains [7, 8] and with BPSV 9108 
(JN162119) [2] and BPSV Iwate/bovine/2007 (AB538385) 
[15]. For detailed phylogenetic analyses, sequences of 12 
previously reported virus isolates from Japan [8] and other 
countries [2] were used (AB044795-AB044801, JN171854, 
JN171860, JN171861, JN191576 and AY386265). 

Virus isolation was performed from 10 swab samples 
as described above. Each sample was mixed with Eagle's 
minimal essential media (MEM) and then centrifuged for 30 
min at 2,500 x g. The supernatant was inoculated onto bo- 
vine testis (BT) cells. The inoculated cultures were observed 
daily for at least seven days after inoculation. Cultures with 
no observable cytopathic effect (CPE) were passaged twice 
more in a blinded manner. After the virus was isolated, we 
performed PCR amplification of the envelope gene followed 
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Drd I Xmn I PflM\ Hinc II 



MABCDEFGHNM MABCDE FGHNM MABCDE FGHNM MABCDEFGHNM 




Fig. 1 . RFLP analysis of eight amplicons. Restriction endonucleases used were Drd I, Xmn I, PflM I and Hinc II. Amplicons A-E and H were 
digested with Xmn I and PflM I, respectively. Amplicons F and G were digested with Xmn I and Hinc II. Lane M, 1 00 bp DNA ladder; A-H, 
digested PCR products; N, undigested amplicon (594 bp). 



BPSV 1 TGGTGTTGTTGGCGGCGTCGTCGCGGCCGGGGATGGAGAAGAGCCGCGTGGAGATGTCGATGCTCCCCACGCC 


GAAGTCGTTC 


AGCGTGCGCGCGGCCGA 


Iwate/bovine/2007 

1 . 










▲ 




A 1 






R 1 






C 1 














D 1 A 






Him II 








E 1 










■ Xmn I 




F 1 


f*i 


A. . 






G 1 


c 


A 






9108 1 T 


c 


A. . 






H 1 


. . . . C. . T. . CAC. . . CTTG. . . AC. . . C. . 


..c. c 


. . G 


G 


. . AC 


PCPV 1 


C ' CAC. . . CTTG. AC C 


..c.c 


. . G 


G 


. . GC 


PVN7 1 fi 


C 


C C 




G 


A f 


G C 


BPSV 101 CACGGACAGCGGGTCCGCGTTCTTCCACTCCGTCACGATCACGCGGACACGCACGCTGCGGTCGATGGCCGCGCGCAGCAACGCGTCCATGATCCGCGGC 
Iwate/bovine/2007 

101 G A 


A 101 T 




G 










B 101 




G 




.... A 






C 101 T 




G 










D 101 T 




G 










E 101 




G 










F 101 . G 




C. G 


0 




C . 




G 101 . G 


PilMl 


C. . G 


c. ... 




C. . 




9108 101 .G 




C. . G 


c. ... 




c. . 




H 101 G.GC. .G 


gT. gg 


r. . . G. T. . . C. . G 


c. ... 


T 


G. . 


G T 


PCPV 101 G.GC. .G 


G. .GG 


r. . . G. T. . . C. . G 


0 


T 


G. . 


. ...G T 


PVNZ 101 . . . C. . G A 


...T. .G. ... 


C. . G 


c. ..c 


. . . c 


G. . 


GGCA. . . 



100 

100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 
100 



200 

200 
200 
200 
200 
200 
200 
200 
200 
200 
200 
200 
200 



BPSV 



201 CAGTACTCCACGCGGTCAGCGTGCCGGATCACGGGCACCAGGGACAGCAGCGACAGGTCGATGCTGTTCTTGGCGGAGTCGATGCGGTGCAGCACCAGGT 300 



Iwate/bovine/2007 
201 ... 



A 
B 
C 
D 
E 
F 
G 

9108 
H 

PCPV 
PVNZ 



201 
201 
201 
201 
201 
201 
201 
201 
201 
201 
201 



. . . GC. CTG TT C. 

. . . GC. CTG TT C. 

. T. . T. CTG G 



. . . TC. . G. 
. . . TC. . G. 

. G. TC. . . . 



. . . . G. 
. . . . G. 
. A. . G. 



. ...T. 
. . . . C. 
. . . . C. 
. A. . C. 
. G. . C. 
~GAC] . 



. CTC. . . 
. TTC. . . 
A. . . . T. 



300 
300 
300 
300 
300 
300 
300 
300 
300 
300 
300 
300 



BPSV 301 CCTCGTCCAGCGTGCGGTAGAAGCCCAGGAACTTCTCGGGCGCGTCCGAGAAGAACACGCCGCCGCCCGCATGGTCCAGGTGGAAGTCCGTCGCCGTTGG 
Iwate/bov ire/2007 

301 T G 


A 


301 




T 










. . . G. . 


B 


301 




T 










. . . . G. . 


C 


301 




T 










fi 


D 


301 














. . . G. . 


E 


301 








G. . 




T. 


. . . 6. . 


F 


301 




G 


G 


G 


. .T 




...G. . 


G 


301 




G 


G 


G. . 


.T 




. . . G. . 


9108 


301 




A. . . . 


. ...G 


G. . 


. .T 




. . . G. . 


H 


301 


. . G 


GCG. . . C A 


. ...G 


G. AG. . 


.T 


G. 


...G. . 


PCPV 


301 


. . G 


GCG. . . C A 


. ...G 


G. AG. . 


.T 


G. 


. . . G. . 


PVNZ 


301 


A 


T 


. ...G.T. .. 


...T G. . 


. . TG. . . . 


A G G 


G. . G. . 


Fig. 2. 



















400 

400 
400 
400 
400 
400 
400 
400 
400 
400 
400 
400 
400 



MOLECULAR ANALYSES OF PARAPOXVIRUSES 1401 



Table 1 . Percent nucleotide and deduced amino acid sequence identities for partial regions of the viral envelope 
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A to H, sequenced amplicons in this study; BPSV, bovine papular stomatitis virus strain V660; PCPV, pseudocowpox virus strain VR634; PVNZ, 
parapoxvirus of red deer in New Zealand strain DPV 



by RFLP as described above. 

Parapoxvirus DNAwas detected from eight swab samples 
(A— H) from six affected calves (farms 1—4, 6 and 8) and 
two calves without clinical symptoms (farms 5 and 7). Five 
amplicons (A-E) were cut with Xmn I only (496 and 98 bp 
fragments), and amplicon H was cut with PflM I (443 and 
151 bp). Amplicons F and G were digested by both Xmn I 
(496 and 98 bp) wdHinc II (516 and 78 bp) (Fig. 1). 

Nucleotide sequences of the PCR products and parapox- 



virus reference strains are shown in Table 1 . Two amplicons 
(E and G) were derived from two calves without clinical 
symptoms. Because the nucleotide sequences of these two 
amplicons were different from those of the other 6 amplicons 
that were derived from calves with clinical symptoms, we 
ruled out the possibility of contamination. 

The location of the Xmn I site was conserved in BPSV 
V660, Iwate/bovine/2007, 7 of the amplicons identified in 
this study (A-F) and BPSV 9108 (Fig. 2). The /fine II site in 



BPSV 401 CGTGATCATGGCGCAGCACAGCCGCGAGAAGGGGACCTTGGGCTCCACGATCATGGAGTAGAAGGTGTTGTACCGGTTCATCAGGTCCGTGGCCAGGTGC 500 
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Fig. 2. (Left page and this page) Nucleotide sequence alignment of partial regions of the envelope gene. Nucleotides identical to BPSV 
are represented by dots. A-H, sequenced amplicons; BPSV, bovine papular stomatitis virus strain V660; Iwate/bovine/2007, BPSV 
strain Iwate/bovine/2007; 9108, BPSV strain 9108 strain; PCPV, pseudocowpox virus strain VR634; PVNZ, parapoxvirus of red deer 
in New Zealand strain DPV. 
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F and G was located upstream of the Xmn I site, because of 
the T61C nucleotide substitution. This substitution was also 
present in BPSV 9108. In the deduced amino acid sequences 
for F and G, the presence of an 130 IV substitution has the 
potential to confer a change in antigenicity. Generally, the 
virus envelope antigen is highly variable, because of its ex- 
posure to host immunological pressures. These substitutions 
might be the result of exposure to some as yet unidentified 
immunological pressures. The nucleotide sequence of ampli- 
con B was identical to that of BPSV Iwate/bovine/2007 de- 
tected in a calf from Iwate Prefecture in 2007 [15]. Based on 
the phylogenetic tree we constructed, seven of the detected 
amplicons (A-G) could be considered BPSV; however, the 
sequence of amplicon H was most similar to that of PCPV 
(Fig. 3). Within the larger cluster, F and G formed a subgroup 
along with BPSV 9108. The nucleotide and deduced amino 
acid sequence identities for the amplicons and parapoxvirus 
reference sequences are shown in Table 1 . 

Three swabs from three calves (farms 1, 2 and 6) showed 
CPE in BT cell cultures after a second blind passage. Other 
swabs showed no signs of CPE. The molecular characteris- 
tics of three isolates were identical to those of amplicons (A, 
B and F) from the swab of the same calves. 

Parapoxvirus envelope genes were detected in eight 
calves from eight farms. Six of these calves were symptom- 
atic, and two were asymptomatic. Because the asymptomatic 
calves had been raised with one symptomatic calf at each 
farm, the calves might have been subclinically infected with 
the virus. Parapoxvirus infections with clinical symptoms 
in cattle have rarely been reported in this prefecture [15], 
although a high rate of antibodies positive against parapox- 
virus has been found [14]. The increase in the level of 
detection is likely a consequence of increased surveillance 
by veterinarians and farmers following a foot-and-mouth 
disease outbreak in the spring of 2010 in Japan. This would 
suggest that parapoxvirus infections are widespread among 
cattle in Iwate. Infections occur irrespective of age with 
relapses common, because infections confer no significant 
immunity [5]. However, clinical symptoms were only ob- 
served in calves from 2- to 10-month-old. This finding is in 
accordance with a previous report that found that infections 
were more common in calves than adult cattle [5]. 

Five (A-E) of the eight samples examined in our study 
were classified as BPSV, based on amino acid identities with 
the corresponding regions in BPSV V660 [8]. In particular, 
amplicon B showed 100% nucleotide identity to the ampli- 
con from Iwate/bovine/2007, in spite of being isolated in 
different years. This strongly suggests that this examined 
region encoding an envelope protein is highly conserved. 
By contrast, the presence of BPSV, which can be digested 
with Xmn I and Hinc II, was also indicated in this study. It 
was reported that BPSV 9108, isolated from a calf in France, 
contained a greater number of nucleotide substitutions in the 
envelope genes compared with other BPSV species [2]. Our 
phylogenetic analysis indicated that the French isolate could 
be grouped together with samples F and G (Fig. 3) and con- 
tained the same substitution at nucleotide 61 (Fig. 2). Our 
results also suggest that BPSV demonstrates a certain level 
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Fig. 3. Phylogenetic tree derived from the deduced amino acid 
sequences of partial envelope gene regions. A-H, sequenced am- 
plicons; V660, BPSV strain V660; VR634, pseudocowpox virus 
strain VR634; DPV, parapoxvirus of red deer in New Zealand 
strain DPV; others, isolates from Japan and other countries. 



of variability. Although there is no direct epidemiological 
evidence, these results suggest that BPSV variants are likely 
widespread among cattle in Japan. 

Sample H was grouped with PCPV, based on RFLP re- 
sults, nucleotide and amino acid identities and phylogenetic 
analysis. This is the first report describing PCPV in Japan. It 
was detected from the lesion around the mouth, not on the 
teats. The present result indicates that disease names are not 
always coincident with causal viral names in parapoxvirus 
infections as previously reported [8-10, 17]. PCPV VR634 
was isolated from a human in the U.S.A. that had "milker's 
nodules" on the hands [4]. The high level of identity between 
sample H and PCPV VR634, whose origins are significantly 
different geographically, biologically and chronologically, 
is indicative of high levels of conservation in the envelope 
gene region we examined. 

Recently, several case reports have highlighted the pres- 
ence and wide distribution of parapoxvirus infections in 
cattle [9-11, 16]. Further molecular studies of parapoxvi- 
ruses circulating among field animals are needed to clarify 
their epidemiology and pathogenesis. 
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